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Introduction 


Assumptions of the model 


Assumptions of the model 


• The irradiation of biological systems leads to the formation of 
radiolytic species such as H-, OH, H 2 , H 2 0 2 , e _ aq , etc. [1] 

• These species react with neighboring molecules, which 
result in damage in biological molecules such as DNA. 

• Radiation chemistry is therefore very important to 
understand the radiobiological consequences of radiation [2]. 

• In this work, we discuss an approach based on the exact 
Green Functions for diffusion-influenced reactions which may 
be used to simulate radiation chemistry and eventually 
extended to study more complex systems, including DNA. 

Green functions 


• The pair of particle may react as follow [3]: 
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k i> k e : product formation rate constant 

• Transitions from a state to another are defined: 
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• The exact Green functions for an isolated pair are 
known analytically [3-4]: 


Sampling of the Green functions 
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• We have developed exact algorithms to sample the 
random variates r for p(r,t|r 0 ) and p(r,t|*) [4]. 

• The algorithm allow the simulation to be done in 
several timesteos (time discretization) 


p(r,t|*): 


k d 

a 

4n:rRD 

(P - a )(Y - a ) 


W 


r-R 

V4Dt 


,-aVDt 


+ ■ 


P 


(y - P)(a - P) 


W 


r-R 

V4Dt 


,-pVBt 


+ ■ 


y 


(a - yxp - y) 


w 


r-R 


,-yVbt 

V4Dt j 


Q(t I *) = 


. k d 

1 (Ra + 1) 

RD 

i ^ 

P 

AO 

-X 

H 

P^ 

AO 

1 

£ 

1 

£ 




P (**,t | *) = 


a(|3 + y) 


(P-a)(y-a) 


Q 


— a-\/Dt)- 


P(y - P)(«. - P) 
P(a + y) 


(y-(3)(a-P) 


n 


-pVot)- 


y(a - y)(P - y) 

y(a + P) 


(a-yXP-y) 


n 


-yVDt) 


k e 

a + P + y 

(P + Y) 

D 

aPy 

P^ 

AD 

1 

P^ 

O' 

1 

P^ 


n(-aVot) n — n(-pVDt) la + |i| — nf-yVot 

v ' P(y-P)(a-P) v ' y(a-yXP-y) V 


The coefficients a, [3 and y are 
rate constants (k D =47iRD): 


related to the reaction 
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W(x, y) = exp(2xy + y )Erfc(x + y) 
f2(x) = exp(x 2 )Erfc(x) 



Left: Green function 47rr 2 p(r,t|r 0 ) for R=1, r 0 =1.5, for k a =47rRD, k d =1 and k e =1 at t=1, 2, 4, 
8, 16 and 32. Analytical functions: ( — ); Result of sampling: (■). Right Survival probability 
Q(t|r 0 ), binding probability p(*,t|r 0 ) and reaction probability p(**,t|r 0 ) as function of time for 
R=1, r 0 =1 .5, k a =47iRD, k d =1 and k e =0. 1 ( — ), k e =1 ( — ) and k e =10 (•••). 


The pair of particles is initially at a distance r 0 and at they 
are at the distance r at t: 

4rcr 0 2 p(r, 1 1 r 0 ) = 8(r - r 0 ) 

The distance between particle obeys a diffusion equation: 
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The material balance condition is: 
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The boundary condition is 
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Survival and dissociation probabilities 
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Many-particles system 

• When more particles are added to the system, the 
number of interactions grow quickly 

• 2 Particles 
-1-2 (1 interaction) 

• 3 Particles 

-1-2, 1-3, 2-3 (3 interactions) 

• 4 Particles 
-1-2, 1-3, 1-4, 2-3, 2-4, 3-4 (6 interactions) 

• N Particles 

-N(N-1)/2 interactions ^ Grows as ~N 2 ! 

• The Green Functions can be used to build a radiation 
chemistry code [4], by using average positions generated 
by sampling the inter-particle distance at each timestep 



Chemical reactions and radiolytic yields 

• The chemical reactions between radiolytic species 
with no electrostatic interaction (i.e. their charge 
product is 0) can be simulated by using the Green 
Functions described above. 

• The radiation chemistry code can be used to 
simulate the time evolution of the radiolytic species 
(radiation chemistry) and radiochemical yields [5,6]. 
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Reaction rate constants (k obs , k dif and k act ), reaction radii (R), probability of geminate 
recombination, and a for reactions between radiolytic species [5]. 


Radiation track structure and evolution in time 





Time-dependent yields of chemical species produced by 300 Time evolution, in 3D, of a 24-MeV 4 He 2+ , LET-26 keV/pm, 

MeV protons (LET -0.3 keV/pm) [6]. Calculations: IONLYS- a t 10' 13 , 10" 9 , 10~ 7 and 10" 6 s. Each dot is a radiolytic species 
IRT ( — ); SBS ( ); The dots are experimental data. 


Conclusion 

• This approach has been used successfully to 
simulate the time evolution of radiolytic species and 
to calculate radiochemical yields. 

• The radiation track structure code RITRACKS [7] 

and the chemistry code will be of crucial 

importance in future models of DNA damage. 
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